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q . In General Relativity, there is a new field of activity concerning the study of charged stars. 

S-h ' In a recent paper, Ray et al. have shown the possibility that the collapse of a charged star 

, could form a charged black hole before all the charge leaves the system. In this field of view 

we propose a new model for UHECR and we will show that it is possible to accelerate cosmic 
rays up to EeV. In this talk we will compute the UHECR flux, the charged black hole density 
and the energy spectrum associated with them in order to reproduce experimental data. We 
■ will see that we need only a small number of these hidden sources in order to account for 

observed UHECR properties and we will study the limits on charge and mass of black holes. 

1 The charge on astrophysical objects 

As it was known since 1924, every self-gravitating system has a net charge. This can be obtained 
moving from the classical Poisson's equations^. For isothermal equilibrium, densities can be 
expressed in terms of Boltzmann distributions, which yield, in terms of expected values for 
energy ^ — etfj ~ Am p (p + Zetp => etp ~ —Am p (j)/ {Z + 1), or in terms of charge to mass density 
ratios p e /p = GAm p /(Z + l)e. For a classical gravitational potential 4> = GM/R and electric 
potential -0 = Q/R- So the net charge of a star of mass M is about 150 C. 

The potential differential between the center and the edge of a self-gravitating system is^ 

A m„ GM A MR™ 



For a neutron star this would be V ~ 10 8 Volts. We notice that a self-gravitating object built 
with particles of extremely high A/Z ratio would have a huge charge. 



2 General Relativity Limits: compact stars and black holes 



Once understood that any self-gravitating object has a charge, here we want to evaluate the 
upper limits on charge that arise from general relativity. 

In radial coordinates, we can write the metric of a spherically symmetric ball of charged 
perfect fluid ds 2 = —e u ^ r ^dt 2 + e x ^ r '^dr 2 + r 2 (d6 2 + sir? 9dcj) 2 ). The stress energy tensor is 
T^ v = (p+P)u^u v + Pg^ v + \tt \F^ a F v a - \g^ v F a P F a& , with 4-velocity, g» v metric tensor, 
p mass density, P pressure and F^ v electro-magnetic tensor. Solving the Einstein equation 
together with the spherical symmetry condition and the Maxwell equations, one obtains the 
hydrostatic equations for a compact charged star, and notably a generalized Oppenheimer- 
Volkoff equation^ 

dP q dq Airr^P + mr — q 2 

dr 47rr 4 dr r(r 2 — 2mr + q 2 ) 
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Coulomb term Regeneration term 

where q(r, t) = Jq e^ x+u ^ 2 4Trr 2 j°dr is the charge inside a shell, and j° is the zero component of 
current density j^. The charge is thus playing a double role, because it contrasts the gravitational 
attraction as well as it contributes to the increase in gravitational mass. So it arises a constraint 
on upper limits on charge. 

The complete structure of a spherically symmetric, charged, compact star was calculated for 
a polytropic equation of state and a charge proportional to density^. It gives Q ~ 10 20 M/Mq C. 
The same equation was solved for a stable quark star, with an EOS p = 3P + 4B, that is the so 
called MIT bag model, and it was found a huge surface field^-^ £ = 10 19 V/m. 

Other constraints can be fixed on black holes. A charged static black hole is given by 
Reissner-Nordstrom metric, that defines the horizon radius as r± = M ± \J M 2 — Q 2 To avoid 
naked singularity, we need to stay below 

Q M = V^GM ~ 1.7 x 10 20 -^- C (3) 

Mq 

A large electric field can be induced by an external magnetic fieldlS] B ~ E 
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J C. (4) 



A maximally rotating black hole, described by a Kerr metric, in a magnetic field aligned 
along its symmetry axis will accrete charge until it becomes^ 

Charge was treated as a perturbation of Kerr metric and cannot force the Reissner-Nordstrom 
limit. Other calculatio ns were made on this, with different magnetosphere structure, finding 
similar values of charge^l 

QED limits: Pai r cr eation puts a strong limit on electric field at surface of compact stars. 
The critical field is^^£ c = m 2 c 3 /?ie = 1.2 x 10 18 V/m and the relative magnetic field limit is 
B c = 4 x 10 13 G. If £ = Q/Qm-, for a black hole the upper limit on electric field is 

S(r+) <£ c ^-^->6- 10 5 - j =- (6) 

1 ; m (i + ^~e) 2 v ; 



Table 1: Potential differential at surface for different stellar objects. Maxima are corrected for QED limits 



Type Min Max 

Normal star 1900 V 7 x 10 2ti V 

Compact star 10 8 V 2 x 10 22 V 

Black hole 10 9 V 1.5 x lO 21 ^ V 



3 Charged black holes as UHECR accelerators: energy spectrum derivation, time 
duration and flux 

If an extremely charged black hole exists, for instance soon after a collapse, one of its signature 
can be the UHECR spectrum. Here we propose a simple model to evaluate the emission of such 
objects. A charged black hole is surrounded by an Hydrogen cloud, which is the reservoir for 
particle to be accelerated. The cloud is ionized if the electric field exceeds the ionization field 
for a H atom: £ > Vi on /ao, Vi on = 13.6 V, ao = 0.5 x 10~ 10 m. The acceleration zone is then 
constrained between the horizon and an outer radius, where the electric field falls below the 
above limit. The extracted protons are accelerated at energy comprised between the following 
limits: 



E max ~ -r^— ^5£ x 10 26 eV ; E mm ~ M^- 6.4 x 10 20 eV (7) 
47re r+ y M 

To calculate the spectrum, we made three simple assumptions: the matter surrounding the 
black hole is the reservoir for UHECR; each remnant atom has the same ionization probability; 
spherical symmetry of remnant. The particle number contained in a radius R is N(R) = 

4:7rr 2 n(r)dr where n is the remnant number density. The particle energy gain is proportional 
to the scaling of potential, E oc r _1 . So the spectrum is 

^-dE = 4Trn(r)( 3 E- 4 dE (8) 
dE 

wh ere C accounts for the efficiency of this process. We know that a self-gravitating acretion disk 
has^^ n(r) oc r 1 ^ 7-1 ) where 7 is the usual polytropic index. Then we can write for spectrum: 

dN . 5-4 7 

—=dE = AE 7-1 dE (9) 
dE 

For a monoatomic gas, 7 ~ 5/3 ==> a ~ 2.5 

Also the duration of emission can be calculated as a toy model, under the assumptions of 
(i) matter density constant uq inside radius ro, (ii) stationary supplied matter in time r and 
(in) constant mass during discharge. So charge will decrease of 

Q = (10) 

That is f = -£ 3 / 2 /T, with Q = £Q M and 

T - 1 = 4JL^L(M_\ 1/2 = M . ltr5 i (M_y* 

3 N t \M C J t \M 3 J 1 ' 

The characteristic dimensions are then = 10 39 , the charged particle number rj on = 2.5T0 9 m, 
the outer ionization radius no = 10 6 m 3 , that is the usual ISM density. The solution is 

£ = £o(l + |^) . (12) 
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Figure 1: Temporal evolution of charge £ = Q/Qm in a black hole and spectrum of an homogeneous distribution 

of such sources, calculated as a MonteCarlo propagation. 



So the discharge time is At = 3 • 10 8 r (M/M ) _1//2 s. We expect r > 10 s (time requested for 
electron to fall inside horizon after ionization), so UHECR bursts duration is At ~ 1 -r- 1000 yr 
The flux at 10 20 eV is Q b s = 3 • 10 -16 m -2 sr~ x s _1 . The flux incoming from a single source 
at Rio = 10 Mpc is 



* = N 2 = 2.8 x 10-"^- 



'm 



sr 



(13) 



10 3 in a GZK volume can account for 
~ 10~ 5 Mpc" 3 (within 95% CL of AGASA 



Few black holes with M/M Q = 10 2 and Q/Qm 
present UHECR flux, that is a source density 

source density^D). Expected active sources may be obtained for instance from local GRB rate 
(~ 10~ 8 Mpc -3 yr -1 for isotropic emission) integrated through lifetime. This result does not 
depend on the process efficiency £, because we exactly know how many particles have to be 
accelerated during the time At. 

As a conclusion, we can say that a charged black hole can accelerate protons up to extremely 
high energy, with a power spectrum according to present observations and well focused on high 
energy. The observed multiplets can also be explained through the long emission time. Possible 
further traces of such electric black holes could be the gamma ray bursts, whose rate is able to 
recover the observed flux. 
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